Tomographic reconstruction of quantum states in spatial dimensions 
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Most quantum tomographic methods can only be used for one-dimensional problems. We show 
how to infer the quantum state of a non-relativistic A'^-dimensional harmonic oscillator system by 
simple inverse Radon transforms. The procedure is equally applicable to finding the joint quantum 
state of several distinguishable particles in different harmonic oscillator potentials. A requirement 
of the procedure is that the angular frequencies of the A*' harmonic potentials are incommensurable. 
We discuss what kind of information can be found if the requirement of incommensurability is not 
fulfilled and also under what conditions the state can be reconstructed from finite time measure- 
ments. As a further example of quantum state reconstruction in A'^ dimensions we consider the two 
related cases of an A^-dimensional free particle with periodic boundary conditions and a particle in 
an Ai'-dimensional box, where we find a similar condition of incommensurability and finite recurrence 
time for the one-dimensional system. 

PACS numbers: 03.65.Wj 



I. INTRODUCTION 



"What do we know about the state of a physical system 
given a certain set of measurements?" . While this ques- 
tion permeates all of physics it is particularly slippery in 
quantum physics. Here, measurements disturb a system 
and potentially alter the outcome of subsequent measure- 
ments. This obstacle is overcome by using a vast ensem- 
ble of uncorrelated and identical quantum systems, where 
one measurement is performed on each ensemble member 
whereafter this member is discarded. The ensemble must 
be chosen large enough to permit measurements of all the 
quantities of interest and to obtain statistically signifi- 
cant data about these quantities. Introducing ensembles 
this way furthermore cements the role of the quantum 
state as the complete statistical information of the quan- 
tum system. Following this ensemble approach, we shall 
in the present paper consider the quantum state as being 
described by the density operator p, which can describe 
both pure and mixed states. Our goal shall be to find 
this operator. 

In principle, since the quantum state is completely 
characterized by its density operator's matrix elements 
in a complete basis {|A)}, one could just measure (the 
real and imaginary values of) all these matrix elements 
(A|p|A') - i.e. the density matrix. However, such a gen- 
eral set of observations may be very difficult to perform 
experimentally, and instead it will be our aim to find the 
quantum state from experimentally realizable measure- 
ments. Specifically, we shall adopt the quantum tomo- 
graphic approach where only measurements of the spa- 
tial distribution is made at different points of time 
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From these diagonal elements of the density operator in 
the position representation, Pr(a;,i) = t{x\p\x)^, and the 
known time evolution due to the Hamiltonian we obtain 
the full density matrix {x'\p\x) or, equivalently, the phase 
space distribution W{x,p). 

Sofar, most methods in quantum state tomography 
have been concerned with systems with only one spa- 
tial dimension. We shall in this paper present the ex- 
tension of two well-known methods of quantum state to- 
mography to N dimensions. In section ^ we shall con- 
sider the harmonic oscillator by a treatment similar to 
that in 0. This is not a trivial extension as revealed 
by a simple consideration of the dimensionality of the 
sets of measurements and the quantum state: Pr(a;, t) 
and {x'\p\x) are both of dimensionality two in the spa- 
tial one-dimensional case, whereas in the iV-dimensional 
case Pr(a:i, . . . , xjq , t) is of dimension TV -I- 1, but the den- 
sity matrix {x[, . . . , x'ff\p\xi, . . . ,xn) is of dimensionality 
2N. In section HTll we will treat the case of free parti- 
cles considered in one dimension in , but with periodic 
boundary conditions and in a box potential. We finally 
give a summary of the paper in section Hvl 



II. THE HARMONIC OSCILLATOR 

It is shown in 6] that there is a 1 : 1 correspondence 
between the density operator and the quantum charac- 
teristic function W{S,), where ^ is a complex variable [T^. 
This means that instead of directly finding p we may just 
as well find the quantum characteristic function. This is 
traditionally the main trick used in the quantum state 
reconstruction of the harmonic oscillator. Before pro- 
ceeding to the multidimensional case we will briefly reca- 
pitulate this procedure in one dimension. The quantum 
characteristic function can be found from p by: 
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Under the harmonic oscillator Hamiltonian H = 
hijj{a)a + 1/2) with angular frequency uo the ladder op- 
erators in the Heisenberg picture evolve according to 



Letting — ut, the posi- 



tion operator x evolves according to: 
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Please note that we use dimensionless coordinates x and 
p We now make a change of variables in from 

the complex number ^ to two real plane-polar coordinates 
(77,0) using ^ = 75?7e'^: 



w{ri,e) 



ie 



= Tr fe'^^e^^p 
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It will prove convenient to choose the variables to be in 
the intervals —00 < 77 < 00 and < < tt. Note that the 
position distribution Pr(a;, 9) at the time 6*/^ is a Fourier 
transform of the quantum characteristic function w{ri,9): 

Pt{x,9) = Tt{p5[x{9) -x]} 



= Tr 
1 

2^ 



00 

d'qw{r],9)e"'^''. 



(4) 



By inverting this Fourier transformation, we can find the 
quantum characteristic function i2;(r/, 9) from the mea- 
sured position distributions: 



w{r),9) 



dxPr{x,9) e*''^. 



(5) 



This is the main equation of quantum tomography in 
one dimension: We can find the quantum state (through 
the quantum characteristic function) by observing the 
position distribution for < 9 < tt corresponding to one 
half period of the oscillator. 

For completeness we write the result in terms of 
the Wigner-function; the complex Fourier transform of 
the characteristic function @, The Wigner function 
is a quasi phase-space distribution whose marginals along 
rotated lines arc the measured position distributions: 



Wix.p) 
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|r/| w(r;, 61) e~'''['°"(^^'''+""(^H 



dr] d9 dx' \7j 



Pr{x\9) e'^i^ 
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The multidimensional oscillator 



We will now proceed to show that it is possible to re- 
construct the joint quantum state of a multidimensional 



harmonic oscillator under certain conditions. 

The Hamiltonian is now H ~ ^jLi ^Wj(ajaj + 1/2) 
and we measure the set of N mutually commuting posi- 
tion operators Xj, j = 1 . . . N. For notational simplicity, 
we arrange these operators in a vector x — [xi, . . . ,xn)- 
The iV-dimensional quantum characteristic function is 
now W{^), where ^ is a vector of N complex variables 
As before, we let £,j = '^Vj^^^^ ^ with —00 < r/j < 00 
and < 9j < tt. This yields the A^-dimensional equi- 
valent of lO: 



= Tr { exp 



N 



i=i 



/OO 
d^xPi{x,e)e''>-- 
-00 
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To gain full knowledge of the function w{r],6) we must 
be able to vary the N variables 9j independently of each 
other on the interval < 9j < tt. This is naturally 
not possible in general, since we can only vary the N 
variables comprising through variation of the one pa- 
rameter t. A way to clearly see this restriction is by 
noticing that the present reconstruction scheme relies on 
a Fourier-transformation, which preserves dimensional- 
ity. While we measure the joint spatial distribution in N 
dimensions for different this equals + 1 dimensional 
measurements, while the quantum state (e.g. the charac- 
teristic function or Wigner function) is a 2A^-dimensional 
object. 

It is important to realize what kind of limitations 
are implied by the inability to vary the 9j^s indepen- 
dently. It is always possible to find the quantum states 
of a single degree of freedom, corresponding to tracing 
out all other degrees of freedom. The limitation comes 
about when trying to find the joint quantum state of 
the A^-dimensional system, in particular the correlations 
between the different degrees of freedom and entangle- 
ment. A simple illustration of this is offered by a two- 
dimensional harmonic oscillator with lui = u;2- Let us 
consider measuring the observable {xi{ojt)x2{ojt)). By 
using we find: 

{xi{ujt)x2{iot)) — {[cos{ujt)xi + sm{Lijt)pi\ X 

[cos{ujt)x2 + sm{ujt)p2] ) 

= cos^(u;t) (2:1X2) + sin^ (ujt) {P1P2) + 

cos(wt) sm^Ljt) {xiP2 + Pi2;2) . (8) 

Here one can see that it is impossible by variation of 
t to find the moments (£1^2) and (^1X2), even though 
one can find their sum. This means that even a simple 
two-dimensional Gaussian state cannot be reconstructed 
if u!i — UJ2- In section lTl Dl we shall give a precise method 
to identify which correlations can be found from a cer- 
tain set of data. Among other results we shall see that 
an A^-dimcnsional Gaussian state can be completely re- 
constructed if no two cjj's are equal. 
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We return now to the problem with ^ : The 9j 's are all 
varied through the one parameter t. The obvious solution 
to this problem is to devise some means to vary the ^^ 's 
independently. There are important situations where this 
is indeed possible, e.g. the case of several entangled light 
fields. In this case a full reconstruction may be done, 
regardless of the values of the luj^s, as can be seen in {TJ. 
The independent variation of 9j can here be achieved 
simply by delaying the measurement on the subsystems 
by introduction of, for example, a variable delay line. 

Another possibility for varying the Oj^s independently 
would be to vary the times for the subsystems inde- 
pendently. For this purpose one might use a method 
closely analogous to the so-called "twin-paradox" from 
special relativity. For instance, imagine two spin-0 par- 
ticles in each their one-dimensional harmonic oscillator. 
One may then leave the one subsystem undisturbed while 
the other is accelerated to a relativistic speed. This sec- 
ond subsystem is then allowed to fly along for a while, 
then accelerated back again towards the first subsystem 
and ultimately brought to rest in its original position. 
The time dilation will hereby delay the second subsys- 
tem compared to the first, effectively giving a means to 
independently vary the elapsed time for the two subsys- 
tems. To avoid direct disturbance of the second oscillator 
due to the acceleration, one should accelerate the system 
perpendicular to its direction of mechanical oscillation. 
In principle, this method of exercising control over the 
time of subsystems can also be used for more than two 
subsystems, but presumably with much increased prac- 
tical complication. We note that the idea of relativistic 
time-displacement of subsystems has been suggested in 
and recently applied to entanglement properties with 
highly non-trivial results 0- 



B. Incommensurable frequencies 

Returning to the general problem of full state recon- 
struction by Q, we shall discuss under what circum- 
stances this is indeed possible. By considering all times 
t > and choosing the w/s mutually incommensurable, 
i.e. their ratios are irrational numbers, we can find a 
unique t to reach any 6 as long as 9j/9k with j ^ k 
is an irrational number. To see this, remember that 
9j = [wjtjTT, with [ ]^ being the modulus function with 
respect to tt. The whole scheme can be pictured as let- 
ting N initially coinciding points move around a circle 
with mutually incommensurable angular frequencies: If 
the points coincide at one angle (which we have chosen to 
be = 0), then no pair will ever again coincide at this an- 
gle. The situation is illustrated for TV = 2 in figure ^ A 
small technical detail in this respect is that since we have 
chosen the intervals of < 0^ < tt and — oo < rjj, Xj < oo, 
then each time a 9j surpasses an integer multiple of tt we 
must let rjj — > — 77^- in {Tj). 

We have found the function w(t7, 6) except on the val- 
ues of where two or more 9j/9k, j 7^ fc is a rational 





the case of a two-dimensional harmonic oscillator. The graph 
to the left shows what values of {61,62) can be obtained for 
uj\ = UJ2- More generally, if lu\/ulI2 ~ 01/02 is an irreducible 
fraction smaller than 1, there will be 02 lines in the (61,62)- 
plane, as shown in the middle graph. Moreover, if 1^1/1^2 
is an irrational number, almost the whole &2)-plane will 
be covered as the measurement time T' — > cx3. The values of 
{61,62) not covered are all the values where 61/62 is a rational 
number. Fortunately, it turns out we do not need this set of 
values to exactly reconstruct the quantum state. 



number. Fortunately, this non-available set of 6 values 
has measure zero, and since w(ri,6) is uniformly con- 
tinuous (and thereby non-singular), the inability to find 
w(t7, 9) on a set of measure zero is of no consequence. 
The uniform continuity of wiv, 0) is a consequence of p 
belonging to the trace class Q. 

The price we pay to gain knowledge of the iV- 
dimensional state as compared to the one-dimensional 
case is that we must measure the joint position distribu- 
tion of all coordinates for all times instead of just half 
of the oscillators' period. In actual applications, where 
infinite measurement times are not available, one would 
presumably use frequencies of the N oscillators whose ra- 
tios are rational numbers and measure for the recurrence 
time of the joint system. The frequencies should then be 
chosen so that the 0-space is sufficiently closely sampled 
for a reliable reconstruction. The exact amount of infor- 
mation obtained in such an experiment will be quantified 
in subsection III Ul 

It should also be noticed that this procedure is equally 
applicable to a quantum system comprised of several 
distinguishable non-interacting particles in separate har- 
monic potentials. 

For completeness we give the formula for reconstruc- 
tion of the A^-dimensional Wigner function: 



W{x,p) = lim 



1 



T' 



00 /'OO 



dt d^T] d"x'x 



— 00 'J —00 



2^(27r)^T' Jq 
J7| (-l)^f=i^^(*)Pr {x',[u:t]^) 



N 



exp ^ i^(-l)^^(*)?7j [xj - x'^ cos{[ujjt]^) ~ pj sa\{[uj jt]^)] 

where Fj{t) = Floor(wjf/7r), and the Floor-function 
rounds downwards to the nearest integer. 
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C. An example of realizing incommensurable 
frequencies 

In the above discussion, it was demonstrated that the 
quantum state of an TV-dimensional harmonic osciUator 
could be exactly reconstructed if the N frequencies were 
incommensurable. It may be noticed that the reconstruc- 
tion of the joint quantum state did not require interac- 
tions between the N degrees of freedom. We will now 
give a brief example of how it is possible, by introducing 
interactions between the oscillators, to reconstruct the 
full quantum state when the all frequencies are identical 
and equal to u. We imagine the N oscillators arranged in 
a line, and introduce nearest-neighbor interaction terms 
in the Hamiltonian. We let k < w be a real coupling 
constant: 

N 

Hint = ^ (a]aj+i + a]_^_^dj^ . (10) 

Arranging now the N annihilation operators in a column 
vector a = (di,a2, • ■ ■ ,0'n)'^, we find from the Heisenberg 
equation of motion: 



d 
di 



d2{t) 



/ K 

K LO K. 









d2{t) 


) 


\ dN{t) j 



where the matrix S) is tri-diagonal and real. Since there 
are no terms containing a!|, we can perform a usual or- 
thogonal diagonalization of this matrix to yield N new 
modes, characterized by new annihilation operators d'j. 
The eigenvalues lo'^ of !3 are well-known from e.g. Hiickel 
molecular orbital theory and solid-state physics: 



2k cos 



N 



In this way, the new ladder operators have simple time 
evolutions d'j{t) = a^(0) exp (— iw^t), and we can again 

find position operators x'At) = l/-\/2[a' (t) 



if we choose k so that the A,- 



s are incommensurable. 



So 
we 



can use the reconstruction method from subsection 111 lj| 
Experimentally, one still has only to measure the i^'s 
since the i^ 's are merely linear combinations of these. 



D. Commensurable frequencies and partial 
information 

We have shown above that a complete tomographic re- 
construction of the state of an A^-dimensional quantum 



oscillator system is possible if all the oscillator angular 
frequencies are mutually incommensurable. This natu- 
rally leads to the question of which aspects of the quan- 
tum state can, and which cannot, be obtained from such 
tomographic measurements if some of the angular fre- 
quencies are commensurable. We will seek to quantify 
this degree of information through the moments of the 
ladder operators dj and aj. For this to be meaningful, 
we must assume that these moments are finite. For con- 
venience we shall be considering the Weyl-ordered (i.e. 
symmetrically ordered) products. Letting r and s be 
iV- vectors with non- negative integer components, these 
moments are generally of the form: 



S{r,s) 



' 3- 



where rj e Nq and Sj e {0, 1, . . . , r^}. 



(12) 



fee {1,2,..., TV}. (11) ©: 



That is, all factors in the product S{r, s) is the sum of all 
symmetric permutations of a number Sj of the operator 
dj and a number {vj — sj) of the operator cij. Here we 
let { }w stand for the Weyl ordering, which is the same 
as the sum of all permutations divided by the number of 
terms. For example: 

{aa^lvK = {a^ajw/ = - {dd^ + d^d) . 

A quantum state is completely characterized if all mo- 
ments of the form in (|12|l are specified. In this way we 
are recasting the question of to what extent the quantum 
state can be reconstructed into the question of how many 
of the moments S{r, s) can be found. 

To keep things transparent, we shall initially consider 
only two oscillators with commensurable frequencies and 
only later generalize to the case of N oscillators. We 
shall need the recurrence time for the system, T = 2tt/u; 
so that 9j = LUjt ~ ajLut = a.j9, making Uj a positive 
integer. Since we are measuring the joint a;-distributions 
at different times, it is natural to consider moments of 
these distributions {[xi{ei)Y^ [x2{02)Y^) ll9|. Recalling 



ri r2 

= ^ £ 5(ri,r2,Sl,S2)e*^["^("^-2.l)+a.(r.-2..)] (^3^ 

Si=0 S2=0 

Since the set of functions {e™^,n € N} is linearly inde- 
pendent on the interval [0;27r[ we can find all the mo- 
ments 5'(ri, r2, si, S2) if there are no two of the expo- 
nential functions in the sum H13|) that have the same 
period in 9. Indeed, to find all moments (and not only 
the symmetric ones) of order (ri,r2), we must know all 
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symmetric moments of this and lower order. A pre- 
cise way to state this is that there must be no recur- 
rences in the following lists, where we keep the r^ 's fixed 
in each list and let the s^ 's assume all possible values: 
(ai [r[ - 2s[] + a2 [r'2 - 2s'2])^, with r'j < rj and s'j < r'^ 

Finding these moments can then be done, for instance, 
by Fourier transformation since the aforementioned ex- 
ponential functions are orthogonal. 

Furthermore, one should notice that reconstructing the 
quantum state through the moments of the ladder opera- 
tors, one in principle needs knowledge only of a small but 
finite interval of the angle 0, and not the whole interval 
[0;7r[. A similar result is found in T^. The fundamental 
assumption that allows for reconstruction from any small 
finite ^-interval is that of finiteness of the moments of the 
ladder operators and position operators in (|13() . The ear- 
lier discussed method of state reconstruction via Fourier 
transformation, 10 , does not suffer from this limitation. 
On the other hand not all these moments need be finite, 
only the ones we use in the reconstruction. 

Finally, one should notice that some of the moments 
can always be found, regardless of the value of ai/a2- 
This trivially includes the moments {xYx^2)j since these 
are directly measured, but also the moments {0102) 
and their complex conjugates for any (ri,r2) since these 
moments evolve with the unique largest numerical fre- 
quency. 

We now proceed to discuss the practical usefulness of 
this approach to reconstruction. It is easy to show that 
it is possible to find all moments with ri < a2 and/or 
r2 < Oil, so one may indeed settle for reconstructing mo- 
ments of only low order, e.g. ri -|-r2 < max(Q!i, Q!2). The 
reason for this is both the difficulty in precisely measur- 
ing higher moments of the joint position distribution and 
that finding higher moments of the ladder operators in 
general requires the ability to measure very rapid vari- 
ations in the joint position distribution (see ifT^ ). In 
addition, it is not necessary to measure a continuum of 
angles if one is only interested in moments up to a certain 
order, but only a number of angles equalling this number, 
which makes the procedure practically feasible. In this 
way one only has to solve a number of equations with an 
equal number of unknowns. 

Generalizing the above results to arbitrary N , we let 
T — 2t: /oj be the recurrence time of the system and 9j = 
ajO. The positive integers aj are arranged in an A'^- 
vector a, and we find that to reconstruct the moments 

up to S{r, s) one needs to measure ^IljLi ^^j^ (^j)^ ^"^^ 
that there can be no recurring numbers in each of the lists 
(once more the 's are fixed in each list while assumes 
all possible values): {a - \r' — 2s'])^, with r'^ < rj and 
s'j 5: r'j. Like before, we can in particular always find all 

moments ^IljLi^j^^; since they are directly measured, 
and also ^Hj!;! ^ ^^"^ complex conjugate for all r. 

It is amusing to note that if all the moments of the 
ladder operators are finite and the angular frequencies are 



incommensurable, it is in principle possible to reconstruct 
the full quantum state from measurements made in a 
small but finite time interval regardless of the dimension 
A'^ - an impossible task if just two angular frequencies are 
commensurable. 

One may also remark that an A^-diniensional Gaussian 
state can always be completely reconstructed if no two 
ujj's are equal: The required measurements are the joint 
position distribution for either 4 points of time or any 
finite continuous interval of time. 

Lastly, we make a short comment on a possible strategy 
for guessing the quantum state in the case of commensu- 
rable frequencies. Even though we do not know the mo- 
ments of the ladder operators individually, we still find 
the sum of two or more - this number determined by the 
how many occurrences there are of a particular number 
in the aforementioned lists. One method of guessing the 
state from an incomplete set of data is the Maximum En- 
tropy Principle, due to Jaynes which has been used 
in several reconstruction schemes ^3-^^. The Maxi- 
mum Entropy principle says that in case one has a set 
of data which could have come about due to several dif- 
ferent quantum states, one should choose the state with 
the largest entropy. In our present scenario, this means 
that if we only know the sum of, say, rir^s different sym- 
metric moments, the Maximum Entropy principle would 
ascribe equal Lagrange multipliers to each observable in 
the Maximum Entropy density operator. 



III. THE FREE PARTICLE 

The tomographic reconstruction of the completely free 
particle has been considered in 5], but as stated herein, 
cannot be used beyond the one-dimensional case. Here, 
we instead study the semi-continuous case of the free par- 
ticle with different boundary conditions. In the first two 
subsections IIII Al and IIII 51 we shall study the free parti- 
cle with periodic boundary conditions, also valid for the 
planar rotor . In the last subsection IIII CI we shall 
give a brief treatment of the particle in a box where it 
will be seen that the different boundary condition has an 
important effect on the available information. 



A. The one-dimensional case 

As in the case of the oscillator, we shall first treat 
the one-dimensional case and later extend this to A^ di- 
mensions. In the present case the reconstruction of the 
quantum state will be done through finding the matrix 
elements of p by simple inversion of Fourier transforms 
used in Let us consider a particle with mass m on 
the spatial interval < x < L. The Hamiltonian is 
H = /2m, yielding the eigenstates \n) with energy 
E{n) — hiln^. Here ft — irh/mL'^ and h is Planck's 
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constant. The eigenstate \n) in the x-representation is 
, , , 1 



(X n) 



VI 



^■Kinx / L ^ — iiln t 



(14) 



We can use this to find the position distribution at any 
time: 

Piix,t) = t{x\p\x)^ 



n,n — — oo 



/\^27ri{'n—n)x/L—iQ,{n^ —n'^)t 



(15) 



n,n' — — oo 



It will now be convenient to change variables from n and 
n' to fi = n + n' , An = n — n' . Note that n and An are 
both either even or odd. Changing summation variables 
in this manner yields: 



11 E E + E E 

n— — oo An— — oo n— — oo An— — c 
Gvon oven odd odd 

' n + An n — An\ 

P 1 : , z I X 



oo oo 



Pr(a;,i) 



2-^1 Ajix I L —iQ.7ilS.nt 

e ' e 



(16) 



Our task is to invert this equation to find the matrix 
elements of p. Fortunately, this is quite easy. Notice that 
there are two exponential functions in H16() and both can 
be used with a Fourier transformation to pick out certain 
values of n and An: The first exponential function in x 
contains only n and once n is fixed, the other exponential 
function can be used to select An. 

We let Nt € N and 2T = 2T:/n be the minimum re- 
quired measurement time, which corresponds to twice the 
time a classical particle with the lowest non-zero energy 
(n = 1) would take to traverse the length L. Also, we 
shall choose /3 7^ and obtain: 

'dx e-2"/9-/i / dt e'"^^*Pr(x, t) = 

2NtT J^NtT 

-"^r-i n— — 00 ^ 
parity as P 

n— — 00 
parity as /3 



p 



2 ' 2 



13 and v of same parity. 



(17) 



We have found all elements of the density matrix except 
those for which /3 = 0, i.e. the diagonal in the (momen- 
tum) n-representation. If we were to choose /3 = 0, we 
would always obtain the result of unity in H17|) . as this 



is the same as taking the trace of p in the x-basis. This 
inability to find the diagonal was also pointed out for fi- 
nite observation times for the completely free particle in 
^ , and will unfortunately carry over to the multidimen- 
sional case. This limitation arises because all probability 



densities of the eigenstates are identical: 



= 1/L 



whereby, for example, a thermal state and any pure eigen- 
state |n)(n| have the same probability distribution at all 
times. Even though finding the diagonal of a density ma- 
trix is impossible given the rest of the matrix, one can 
use the Schwartz inequality to constrain the size of the 
diagonal elements through |p(n,n')p < p(n,n)p(n' ,n'). 



B. A'^-dimensional free particle with periodic 
boundary conditions 

We will now move on to the iV-dimensional case. We 



shall work in the TV-dimensional interval Xj € [0,L 



The energies of the eigenstates are E{n) = h ^2^=1 ^. 
and in extension of H16|l we introduce the iV-vectors of 
intergers n and An. This yields: 



2 



Pr(a;,t) 



N 



n EE+EE 

j=l \ rij Arij rij Arij 
even even odd odd' 



^2111 An j Xj j Lj ^ — iflj rij Artjt 



n + An n — An 



(18) 



We shall try to invert this equation to find the matrix 
elements of p. Of course, we cannot approach this com- 
pletely as in the one-dimensional case and select a par- 
ticular vector n through Fourier-transforms, having only 
the parameter t to vary. We can, however, expand the 
time interval of integration to incorporate all points of 
time: 



1 



T" 



T"-»oo 2T" j_rp„ 



(19) 



If the f2j 's are mutually incommensurable, the only possi- 
bility for this Kronecker delta-function to give a non-zero 
result is for = fijij. Remembering that ^Ij = '^'l'2 , 

rrij J 

the condition of incommensurability of the fij's is equi- 
valent to demanding the elements of the list: 



(irijL'j^ be incommensurable 



In this way it becomes possible for a single delta-function 
to effectively serve as N distinct delta-functions. We 
shall henceforth make this requirement of incommensu- 
rability of the rjo 's and, as in the one-dimensional case. 
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we only obtain useful results for f3j 0. With these limi- 
tations we can exactly reconstruct the rest of the density 
matrix. Additionally, Pj and fij have the same parity for 
all j: 



u + f3 u-p 



X X 



dt d'^C 
-T" Jo 



(20) 



The special feature of this system, which allows the treat- 
ment above, is that all its eigenenergies are rational num- 
bers (actually a integers) times some minimum energy. 
This can in fortunate circumstances, i.e. if the products 
of the energy eigenstates in the position representation 
are reasonably placid, allow state reconstruction equa- 
tions like that in H17|) . These reconstruction equations 
will make use of a finite measurement time, and by us- 
ing the trick in (|19|l . can be used in multiple dimensions. 
Another example of this kind of system is the quantum 
mechanical rotor with fixed angular momentum projec- 
tion, which will be the topic of a forthcoming article . 



C. Particle in a box 

In the last two subsections we found ourselves unable 
to determine the momentum distribution for the free par- 
ticle with periodic boundary conditions. The reason for 
this was that all the energy eigenstates of the system had 
the same spatial distribution. In this subsection we shall 
show how imposing other boundary conditions can com- 
pletely alter this situation. In particular, we shall study 
the particle in a box, where the spatial density is zero 
at the boundaries of the box. As usual, we shall first 
treat the one-dimensional case and later move on to the 
TV-dimensional case. 

We choose x G [0, L] whereby E{n) — hil'-n? , now with 
n e N. We use ^' — — yielding the energy 

eigenstates: 



2 . f TIX\ 

— sm n — 
L \ L ) 



-in't 



These eigenstates all have different spatial distributions, 
and we already suspect that we shall be able to recon- 
struct the momentum distribution. Following section 
nil Bl we find the joint position distribution, similar to 



Pr(x,i) = t{x\p\x)^ 

oo 2 



n,n' — 1 
oo 



cos 



{n~\- n) — 



Making once more the substitution n = n + n' and An 
n — n' we arrive at the box equivalent of (|16|) : 



Pr(x,i) 



1 

I 



n— An— — n+1 ^—1 An— — n+1 
Vcvcn even odd odd 



E 



n - 



An n - A 



(^An—^ — cos 



X 

7ra;\ 

t) 



^ — i^nAnt (21) 



Recalling that the set of functions {cos(fc7ra;/L)}, /c G N, 
is orthogonal on a; G [0, L] we can use the cosine functions 
in H21() and the exponential function in time to select a 
certain term in the sum. Selecting G N and /3 G Z with 
> |/3| and letting T' ~ 2TT/il' we obtain the equivalent 

of dni) M- 



1^ + (3 V- [3 



L I 



1 



2NtT' 



2 /da; cos I 



NtT' 

dte'^'''l^'PT{x,t), i'>\p\. (22) 



Here it is apparent that we can reconstruct the full den- 
sity matrix, including the momentum distribution. Thus, 
the choice of boundary conditions has allowed us to over- 
come the limitation encountered in the case of periodic 
boundary conditions. 

Moving on to the multidimensional case where E{n) = 
fiJ2j=i ^^"^ equivalent of H18|) : 



Pr{x,t) 



\ 



E E ■ 

Uj—O Arij — 
even -nj+1 
even 



oo 

-E E 

rij—l Auj — 
odd 



"3-1 / 

odd / 



1 



cos 



An,— ^ 



cos n 



nx^ 

n + An n — An 



^-iUjrijAnjt I ^ 



(23) 



Selecting the two A^- vectors v and (3 with > \(3j\ and 
demanding incommensurability of the elements in the list 
[rrijLj^, we find the reconstruction formula correspond- 
ing to (EOJ: 



u + 13 u^p 



= lim 

T"^oo 2T" 

N , 

ncos(/3,^ 

j=i ^ 

Pr(a;,t), 



T" i-L 

dt d^X X 
-T" Jo 

i^,>m. (24) 
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As in the one-dimensional particle in a box, it is in the 
multidimensional case possible to completely reconstruct 
the density matrix. 

IV. SUMMARY 

We have shown how to extend two common schemes 
of quantum state tomography from one to N dimensions: 
The harmonic oscillator and the free particle on a fi- 
nite interval. In both cases a complete reconstruction 
required extension of the time interval of observation to 
all times and required incommensurability of the eigen- 
energy differences for the different N spatial dimensions. 

For the harmonic oscillator, we quantified the informa- 
tion that can be reconstructed if some of the N frequen- 



cies are commensurable. This was done through recon- 
structing the moments of the ladder operators, and also 
constitutes a reconstruction method for these moments if 
they are finite. This partial reconstruction only required 
measurements at a finite number of points of time. 

For the free particle on a finite interval we showed 
that all off-diagonal elements of the density matrix in the 
energy-representation can be reconstructed in the case of 
periodic boundary conditions, and that a full reconstruc- 
tion is possible for the box potential 
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